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SUMMARY 


A  method  has  Been  developed  for  evaluating  from  shadow  or  schlieren 
photographs  the  pressure  drag  of  axially  symmetrical  Bodies  at  zero 
angle  of  attack  or  of  two-dimensional  Bodies  producing  detached  shock. 

The  method  consists  in  the  determination  of  the  flow  properties  along 
a  characteristic  line  and  the  application  of  the  momentum  theorem  to  a 
stream  tube  around  the  Body.  The  method  can  Be  applied  to  the  determina¬ 
tion  of  external  drag  of  supersonic  inlets  with  suBsonic  flow  at  the 
entrance  of  the  inlet. 


INTRODUCTION 


The  system  of  characteristics  permits  the  pressure  drag  and  the 
shape  of  the  shock  for  axially  symmetrical  and  two-dimensional  flow 
phenomena  having  everywhere  supersonic  speed  to  Be  determined  directly 
from  the  Boundary  conditions.  In  this  case  if  viscous  effects  are 
neglected,  the  value  of  the  pressure  drag  can  Be  determined  analytically, 
and  no  experimental  measurement  is  necessary.  The  determination  of  the 
shape  of  the  shock  can  Be  useful  in  making  a  comparison  Between  experi¬ 
mental  and  theoretical  results  and  permits  an  evaluation  of  the  displace¬ 
ment  thickness  of  the  wake  and  of  the  Boundary  layer  along  the  body. 

The  characteristic  system  does  not  permit,  however,  direct  determination 
of  the  pressure  drag  when  the  shock  is  detached  from  the  Body,  Because  it 
cannot  Be  applied  to  the  subsonic  part  of  the  flow. 

The  determination  of  the  pressure  drag  of  bodies  producing 
detached  shocks  is  important  in  many  problems;  for  example,  the  Bodies 
of  revolution  used  in  practical  applications  often  have  round  noses  and, 
in  supersonic  inlets,  conditions  exist  in  which  the  flow  at  the  entrance 
of  the  inlet  is  suBsonic.  The  experimental  measurement  of  the  pressure 
drag  for  these  cases  usually  requires  complicated  experimental  equipment; 
therefore,  the  determination  of  drag  from  shadow  or  schlieren  photographs 
of  the  phenomena  can  Be  very  useful  for  practical  applications.  A 
method  is  presented  herein  that  permits  determination  of  the  pressure 
drag  directly  from  shadow  or  schlieren  photographs.  This  system  can  Be 
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especially  useful  in  all  cases  in  which  the  experimental  technique  uses 
the  principle  of  moving  models  (shooting  ranges,  whirling  arms,  etc.) 
nnfi  in  experiments  of  supersonic  inlet b  in  which  also  the  internal  flow 
produces  forces  and  in  which  it  iB,  therefore,  difficult  to  measure  the 
external  drag  independently  of  the  internal  drag.  The  method  gives  the 
value  of  the  pressure  drag  due  to  the  subsonic  and  supersonic  regions  of 
the  flow  and  permits  also  the  determination  of  pressure  and  velocity 
distribution  of  the  supersonic  region,  from  which  the  shape  and  thickness 
of  the  wake  can  be  obtained. 


SYMBOLS 


x,  y  Cartesian  coordinates 

n  normal  to  streamline 

Y  velocity 

y  limiting  velocity 

L 

w  velocity  in  terms  of  the  limiting  velocity  (V/V}) 

M  Mach  numb©!’ 

p  density 

p  pressure 

s  entropy 

mg  mass  flow 


l  = 


sin  6  tan  P  sin  cp 


cos(cp  +  p) 

two-dimensional  phenomena 


for  axially  symmetrical  phenomena;  1  =  0  for 


sin  P  tan  P  sin  cp  ^  axially  symmetrical  phenomena;  m  =  0  for 
cos(q>  -  p) 

two-dimensional  phenomena 

inclination  of  the  velocity  vector  with  respect  to  x-axis 
inclination  of  the  shock  with  respect  to  upstream  velocity 


P 


Mach  angle 
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specific  heat  at  constant  volume 
ratio  of  specific  heats 
gas  constant 


Subscripts : 


refers  to  stagnation  free— stream  conditions 
refers  to  free— stream  properties 
points  of  first  family 
points  of  second  family 

quantities  at  the  points  cal  ciliated  from  A  and  B 
quantities  at  any  point  P  behind  the  shock 


THE  METHOD  FOE  DETERMINING  THE  PEESSURE  DEAG 


In  any  supersonic  phenomenon  in  which  a  detached  shock  exists,  the 
flow  is  subsonic  only  in  a  small  zone  of  the  space  and  the  speed  outside 
of  this  zone,  which  is  usually  limited  to  the  vicinity  of  the  nose  of 
the  body,  becomes  again  supersonic.  The  body  considered  has  finite 
dimensions.  In  the  supersonic  part  of  the  flow  the  characteristic  theory 
can  be  applied  to  determine  the  pressure  drag  and  the  pressure  distribu¬ 
tions.  The  line  HN  of  figure  1  is  the  sonic  line  that  divides  the 
subsonic  from  the  supersonic  region,  and  from  any  point  G  of  the  shock 
in  the  supersonic  region  the  characteristic  line  GD  of  th9  second 
family  that  meets  the  body  at  some  point  D  (see,  for  example,  reference  l) 
can  be  determined. 

If  the  flow  properties  are  known  along  the  line  GD,  the  component 
of  the  momentum  of  the  stream  along  GD  in  the  direction  of  the  free— 
stream  velocity  and  the  resultant  of  the  pressure  on  the  surface  GD  in 
the  same  direction  can  be  evaluated.  The  resultant  of  the  pressure 
along  the  surface  OD  in  the  direction  of  the  x-axis  can  be  obtained 
from  the  momentum  law.  Indeed,  the  difference  of  momentum  between  the 
surfaces  EE1  and  GD  must  be  equal  to  the  resultant  of  pressures  at 
the  surfaces  E*G,  EE1,  EO,  OD,  and  GD.  The  pressure  along  EfG  and  EO 
does  not  produce  any  component  of  force  in  the  direction  of  the  x-axis. 

The  resultant  of  pressure  at  the  surfaces  EE*  and  GD  is  known  (along  EE1 
the  pressure  is  equal  to  the  free— stream  pressure);  therefore,  the 
resultant  of  the  pressure  in  the  direction  of  the  x-axis  along  OD  can 
be  evaluated  from  the  difference  of  momentum  in  the  direction  of  the 
x-axis  between  EE1  and  GD.  When  supersonic  inlets  having  subsonic 
entrance  velocity  are  considered  (fig-  2),  the  streamline  EO  is  not 
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parallel  to  the  stream  direction  and,  therefore,  with  this  system  the 
pressure  drag  along  this  line  is  also  determined. 

In  these  considerations  the  viscous  effects  have  "been  neglected,  hut 
the  viscous  effects  are  small  in  the  front  part  of  the  "body  and  therefore 
can  he  considered  independently.  For  the  hack  part  of  the  body  the 
considerations  made  in  the  "Introduction"  on  the  possibility  of  the 
determination  of  the  displacement  thickness  of  the  boundary  layer  and  of 
the  wake  are  still  valid.  It  should  he  noted  that,  by  consideration  of 
the  change  in  momentum  in  the  y  direction,  the  lift  of  two-dimensional 
bodies  at  an  angle  of  attack  can  be  determined.  From  the  preceding  con¬ 
siderations  the  determination  of  the  pressure  drag  along  EOD  (fig.  l)  is 
reduced  to  the  determination  of  the  flow  properties  along  GD.  The  flow 
properties  along  GD  can  be  determined  from  the  shape  of  the  shock  in 
the  following  way: 

If  EGL  is  the  shape  of  the  shock  determined  from  experimental 
measurements  (shadow  or  schlieren  photographs),  the  inclination  of  the 
shock  along  GL  (fig.  l)  can  be  measured.  In  order  to  obtain  the  inclina¬ 
tion  of  the  shock  with  some  precision  the  coordinates  of  the  Bhock  can 
be  measured  and  an  analytical  expression  can  be  determined  for  the  curve. 

From  the  inclination  of  the  shock  e,  the  magnitude  W  and  the 
direction  cp  of  velocity  and  the  variation  of  entropy  As  across  the 
shock  can  be  calculated  at  any  point  P  of  the  shock,  for  example,  from 
the  equations: 


1 

tan  cpp 


7+1 


M- 


\ 


Mj_2sin2e  —  1 


-  1 


tan  f 


(i) 


tan  6p 
tan(e  -  <p)p 


7-1 
7  +  1 


1  —  Wp2coe2(e  —  <p) 


Wp2sin2(e  -  cp) 


(2) 


Pp  tan  6p 
P1  tan(e  -  cp)p 


(3) 
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2  7 

7  +  1 


sin^e 


P 


00 


pp/pl\7 

Asp  =  cv  loge  —  — 

*1  V p/ 


(5) 


In  the  equations,  the  subscript  P  indicates  quantities  at  the 
point  P  behind  the  shock  and  the  subscript  1  indicates  free-stream 
conditions.  From  two  points  of  the  shock  A  and  B  (fig.  l),  the 
flow  properties  at  the  point  C  can  be  determined  by  means  of  the  charac¬ 
teristic  system  for  rotational  flow  as  is  shown  in  the  following 
discussion. 

From  A  the  tangent  to  the  characteristic  line  of  the  first  family 
and  from  B  the  tangent  to  the  characteristic  line  of  the  second  family 
can  be  drawn  (reference  l) 


=  tan(cp  +  p)A 


> 


(6) 


(Si  '  t“(<P  *  13  >B 

j 


where 


sin^(3 


7-l/l 

2  Vw2 


- 1 


The  intersection  of  the  two  characteristic  lines  determines  the 
point  C  in  the  first  approximation.  The  properties  at  C  in  the  first 
approximation  can  be  determined  from  the  following  expressions  (reference  l) 
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—  w  ■—  -  tsn  p(<Pa  -  ^c)  - 


l  4* 


da  _1_  Bln3p  / 
^  XR  cos(p  +  cp/  ^ 


xc)  =  O  (7) 


%  -w( 

w 


tan  p(q>B  -  <PC) 


XB  ~ 


m 


ds  _1_  sln3p  / 
&a  cos(cp  — 


*c)  =  O  (8) 


where,  for  two-dimensional  phenomena,  l  and  m  are  zero  whereas,  for 
axially  symmetrical  phenomena. 


l  = 


sin  (3  tan  p  sin  cp 
cos  (9  +  p) 


(9) 


m  = 


sin  p  tan  p  sin  <p 
cos(qp  —  p) 


(10) 


In  the  first  approximation  the  quantities  with  no  index  in  equation  (7) 
can  be  assumed  equal  to  the  corresponding  quantities  at  the  point  A 
and  in  equation  (8)  to  the  quantities  at  the  point  B,  and  for  ds/dn  the 
following  expression  can  be  used: 


ds 

dn 


_ ^A  ~  _ 

'  sin  PA  /  ^  sin  Pp 

\XA  XC'cos(p  +  cp)^  +  X^'cos(cp  —  p)B 


(11) 


In  this  approximation  the  entropy  gradient  is  assumed  to  be  the  same 
in  the  zone  between  B  and  C  as  in  the  zone  between  A  and  B.  In 
the  process  of  calculations  it  is  possible  that  the  two  known  points 
fall  on  the  same  streamline.  In  this  case,  the  value  of  ds/dn  for  the 
first  approximation  can  be  assumed  equal  to  the  value  of  ds/dn  at  one 
of  the  two  points.  The  value  of  ds/dn  can  be  obtained  from  the  differ¬ 
ence  in  entropy  at  this  point  and  at  a  nearby  known  point  on  the  same 
characteristic  line. 

After  the  properties  at  the  point  C  in  the  first  approximation  are 
obtained,  a  second  approximation  can  be  obtained  in  the  following  way:  The 
value  of  entropy  at  the  point  C  is  determined  from  the  entropy  of  the 
point  on  the  shock  that  is  in  the  same  streamline  as  the  point  C„  This 
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point  is  determined  from  mass— flow  considerations.  The  flow  at  C  is 
then  determined  again  by  averaging  the  values  of  the  known  terms  in 
equations  {6),  (l),  and  (8)  between  the  values  at  A  and  C  and  at  B 
and  C,  respectively. 

For  axially  symmetrical  phenomena  the  masB  flow  contained  in  the 
stream  tube  bounded  by  the  two  surfaces  of  revolution  having  as  generatrix 
the  streamline  that  passes  at  A  and  C  (mA)  and  at  B  and  C  (nig)  can 
be  determined  from  figures  1  and  3: 
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mA  =  /  2itypY 

VC 


A  _  „  sin  p 


sin(cp  +  p) 


dy 


“b  = 


ft B 


J. 


yc 


sin  p 

— 2itypY  - ; - r  dy 

sin(cp  —  p) 


(12) 


The  density  at  any  point  P  can  be  expressed  as  a  function  of  the 
free— stream  stagnation  density  and  of  the  local  velocity  and  entropy  in 
the  form 


As  1 

pp-p^Ci-V)7-1 


(13) 


and,  therefore,  can  he  calculated  at  A,  B,  and  C.  The  velocity  at  any 
point  along  AC  and  BC  can  he  expressed  in  the  form: 


Therefore, 


Y  =  YA[l  +  Ai(y  -  yA)] 
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where 


OAa) 


Similar  procedures  may  he  used  to  obtain 


£  +  (7  ~  ^a)bi] 


sin  (3 

sin(qp  -I-  p) 


[j  +  Cl(y  -  ?A)] 

sin(cp  +  p)A  L.  1  -1 


where 


7C  JA 


(l^h) 


C 


1  = 


sin  pc  sin(cp  +  p)A 
sin  PA  sin(cp  +  p)q 


7C  “  7A 


(14c) 


If,  therefore,  terms  of  higher  order  are  neglected. 


mA  =  *pava  t\y;  (7 a  ~  7c)j  (7c  +  yA) 


sin(cp  +  p)A 


+ 


g7C  +  lA 

3 


+  Bj_  + 


(15) 


MCA  TN  No.  1808 


9 


Similarly, 


where 


yp  Ig(yc  _  7b)(a2  +  B2  +  C2)] 


(16) 


^2 


■CM- 

Vtb 


-  yB 


Bo  =  I  —  —  1 


PB  J7q-  yB 


> 


C2  = 


sin  Pq  sin(cp  -  p)B 


sin  PB  sin(qp  -  p)c 


-  1 


yc  -  yB 


(17) 


If  a  two-dimensional  phenomenon  is  considered. 


mA  =  -pAVA 


sin(q>  !\)a  (yC  "  7k)  +  V‘  -g—  (ax  +  B1  +  cl) 


(yc  “  *aJ 


(18) 


™B  =  PBVB 


/ 

(yc  ^e) 


8ln  H  frc  -  yB>  -  ♦  %  *  %> 


sin(cp  -  p)B 


(19) 
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When  the  mass  flow  is  determined.,  the  points  C*  and  C**  on  the 
shock  along  the  same  streamline  that  passes  at  C  can  he  determined  for 
axially  symmetrical  phenomena  from  the  expressions 


and,  for  two-dimensional 


mA  =  plVl(yA2  “  7C2)* 
“b  =  pivi(yB2  -  yc2)rt 

phenomena, 

mA  =  P1T1(yA  -  yc) 

“b  =  pivi(yB  ■  7c) 


In  the  first  approximation,  Cf  vill  "be  slightly  different  from  C,f; 
therefore,  a  point  Cf,t  vill  have  to  “be  chosen  between  the  two  points. 

The  variation  of  entropy  at  C,,f  can  be  determined  from  equation  (5), 
and  the  entropy  at  Cf?f  behind  the  shock  is  the  same  as  the  entropy  at  Co 
Nov  at  the  point  C,  the  quantities  As,  cp,  P,  and  ¥  are  knovn  in  the 
first  approximation;  therefore,  a  second  approximation  can  be  determined 
by  assuming  for  the  quantities  without  index  in  equation  (7)  the  corre¬ 
sponding  average  values  between  A  and  C,  and  in  equation  (8)  the  average 
values  between  B  and  C.  In  equation  (7)> 


ds 

dn 


(sA  -  S0)jcos(cp  +  p)A  +  cos(9  -  P)q] 

(za  ~  *0)  (sin  Pc  +  Bin  Pa) 


and  in  equation  (8), 


(sB  -  sG)fcos(cp  -  P)B  +  cos(9  -  P)CJ 


ds 

dn 


(xB  -  xc)  (sin  Pc  +  sin  pB) 
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A  second  approximation  can  be  obtained  also  for  tbe  points  C* 
and  C,f  and,  therefore,  for  the  value  of  As  at  C.  If  necessary, 
a  higher  approximation  can  be  calculated  for  As.  Proceeding  in  a 
similar  way  permits  all  the  characteristic  net  CDL  to  be  determined 
and,  therefore,  the  flow  properties  along  CD  can  be  obtained.  For 
inlets  the  calculation  also  gives  the  value  of  mass  flow  that  goes 
inside  the  inlet,  because  it  gives  the  position  of  the  point  E  that 
limits  the  stream  tube  that  goes  inside  the  inlet. 

A  check  of  the  precision  of  the  calculations  can  be  obtained  from 
the  comparison  of  the  mass  flow  contained  in  the  stream  tube  EE’  with 
the  mass  flow  across  CD.  In  addition,  the  entropy  at  D  must  be  equal 
to  the  entropy  at  0  (fig.  l). 

if  the  tests  are  performed  on  a  fixed  model,  a  pressure  measurement 
at  a  point  in  the  zone  DF  can  also  be  used  for  control  of  the  precision 
of  the  calculations.  This  control  is  useful  when  inlets  are  considered 
because  the  stream  tube  E**E’  (fig.  2)  in  this  case  is  not  known.  When 
the  flow  field  along  CD  has  been  determined,  the  pressure  drag  along  EOD 
can  be  determined  from  the  momentum  equation,  and  from  the  shape  of  the 
shock  GL  and  from  the  characteristic  line  GD,  the  pressure  nln-ng  DF 
and  the  shape  of  the  streamline  DT  can  be  determined.  The  streamline  DT 
gives  the  displacement  thickness  of  the  boundary  layer  and  of  the  wake. 


CONCLUSIONS 


A  method  for  determining  pressure  drag  in  axially  symmetrical  or 
two-dimensional  flow  phenomena  from  shadow  or  schlieren  photographs  has 
been  developed.  The  method  uses  the  characteristic  system  for  rotational 
flow  and  can  be  applied  when  detached  shocks  are  produced  by  the  body. 

The  method  can  be  applied  to  the  determination  of  mass  flow  and 
external  drag  of  supersonic  inlets  with  subsonic  flow  at  the  entrance  of 
the  inlet . 
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